V9.3-4 Surface Integrals 



3. Flux through general surfaces. 



For a general surface, we will use ^^-coordinates. It turns out that here it is simpler 
to calculate the infinitesimal vector dS — ndS directly, rather than calculate n and dS 
separately and multiply them, as we did in the previous section. Below are the two standard 
forms for the equation of a surface, and the corresponding expressions for dS. In the first 
we use z both for the dependent variable and the function which gives its dependence on x 
and y; you can use f(x,y) for the function if you prefer, but that's one more letter to keep 
track of. 



(11a) z = z(x,y), 

(lib) F(x, y, z) = c, 



dS = (—z x i — z y j + k)dxdy (n points "up") 
(choose the right sign); 



dS = ± — — dx dy 



Derivation of formulas for dS. 

Refer to the pictures at the right. The surface S lies over its projection R, 
a region in the xy-plane. We divide up R into infinitesimal rectangles having 
area dx dy and sides parallel to the xy-axes — one of these is shown. Over it 
lies a piece dS of the surface, which is approximately a parallelogram, since its 
sides are approximately parallel. 

The infinitesimal vector dS = n dS we are looking for has 

direction: perpendicular to the surface, in the "up" direction; 
magnitude: the area dS of the infinitesimal parallelogram. 

This shows our infinitesimal vector is the cross-product 

dS = A x B 

where A and B are the two infinitesimal vectors forming adjacent sides of 
the parallelogram. To calculate these vectors, from the definition of the 
partial derivative, we have 

A lies over the vector dx i and has slope f x in the i direction, so A = dx i 
B lies over the vector dy j and has slope f y in the j direction, so B = dy j 
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which is (11a). 

To get (lib) from (11a), , our surface is given by 

(12) F(x,y,z)=c, z = z{x,y) 

where the right-hand equation is the result of solving F(x, y, z) = c for z in terms of the 
independent variables x and y. We differentiate the left-hand equation in (12) with respect 
to the independent variables x and y, using the chain rule and remembering that z = z(x, y): 



F(x,y,z) = c 



dx 

' dx 



F, 
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from which we get 



dz F x . <9z F y 

— = - — , and similarly, — = - — 
dx F z dy F z 



Therefore by (11a), 

dS = 
which is (lib). 



dz dz 
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dx dy 
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Example 3. The portion of the plane 2x — 2y + z = 1 lying in the first octant forms a 
triangle S. Find the flux ofF = .xi+yj+zk through S; take the positive side of S as the 
one where the normal points "up" . 

Solution. Writing the plane in the form z = 1 — 2x + 2y, we get using (11a), 
dS = (2i - 2j + k )dxdy , so 
F • </S = / / (2x-2y + z)dydx 

(2x -2y + (l-2x + 2y)) dy dx , 



y 


1/2 x 




R . 
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where R is the region in the xy-plane over which S lies. (Note that since the integration 
is to be in terms of x and y, we had to express z in terms of x and y for this last step.) 
To see what R is explicitly, the plane intersects the three coordinate axes respectively at 
x = 1/2, y = —1/2, z = 1 . So R is the region pictured; our integral has integrand 1, so its 
value is the area of R, which is 1/8. 

Remark. When we write z = f(x,y) or z = z(x,y), we are agreeing to parametrize 
our surface using x and y as parameters. Thus the flux integral will be reduced to a double 
integral over a region R in the xy-plane, involving only x and y. Therefore you must yet rid 
of z by usiny the relation z = z(x,y) after you have calculated the flux integral using (11a). 
Then determine the region R (the projection of S onto the ccy-planc), and supply the limits 
for the iterated integral over R. 



Example 4. Set up a double integral in the xy-plane which gives the flux of the field 
F = x i + yj + zk through that portion of the ellipsoid 4x 2 + y 2 + Az 2 = 4 lying in the 
first octant; take n in the "up" direction. 

Solution. Using (lib), we have dS = ^—^^^-dxdy. Therefore 

8z 

[[F.dS= [[ 8x2 + 2 f + 8z2 dxdy= ff 1 -dxdy= (I dX "y , 

JJs JJs 8z JJgz JJ R ^l-x 2 -(y/2) 2 

where R is the portion of the ellipse Ax 2 + y 2 = 4 lying the the first quadrant. 

The double integral would be most simply evaluated by making the change of variable 
u = y/2, which would convert it to a double integral over a quarter circle in the xu-plane 
easily evaluated by a change to polar coordinates. 



V9.3-4 SURFACE INTEGRALS 



4. General surface integrals.* The surface integral J J f(x, y, z) dS that we 
introduced at the beginning can be used to calculate things other than flux. 

a) Surface area. We let the function f(x,y,z) = 1 . Then the area of S = J J dS . 

b) Mass, moments, charge. If S is a thin shell of material, of uniform thickness, and 
with density (in gms/unit area) given by 5(x,y,z), then 

(13) mass of S = S(x,y,z)dS, 



(14) x-component of center of mass = x = f [ x ■ S dS 

mass S J Js 

with the y- and z-components of the center of mass defined similarly. If S(x, y, z) represents 
an electric charge density, then the surface integral (13) will give the total charge on S. 

c) Average value. The average value of a function f(x,y,z) over the surface S can be 
calculated by a surface integral: 

(15) average value of / on S = f(x,y,z)dS . 
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Calculating general surface integrals; finding dS. 

To evaluate general surface integrals we need to know dS for the surface. For a sphere 
or cylinder, we can use the methods in section 2 of this chapter. 

Example 5. Find the average distance along the earth of the points in the northern 
hemisphere from the North Pole. (Assume the earth is a sphere of radius a.) 

Solution. — We use (15) and spherical coordinates, choosing the coordi- 
nates so the North Pole is at z — a on the z-axis. The distance of the point 
(a, <p, 9) from (a, 0, 0) is a<p, measured along the great circle, i.e., the longi- 
tude line — see the picture). We want to find the average of this function 
over the upper hemisphere S. Integrating, and using (9), we get 

a<j>dS = I / acf>a 2 sin <f>d(j>d9 = 2na 3 / (f> sin (f> d(j) = 2-rra 
J Js Jo Jo Jo 

(The last integral used integration by parts.) Since the area of S = 2-kcl 2 , we get using (15) 
the striking answer: average distance = a . 

For more general surfaces given in xyz-coordinates, since dS = ndS, the area element 
dS is the magnitude of dS. Using (11a) and (lib), this tells us 

(16a) z = z(x,y), dS = z% + z 2 + l dxdy 

IVFI 

(16b) F(x,y,z) = c, dS=——dxdy 

\r z \ 

Example 6. The area of the piece 5* of z = xy lying over the unit circle R in the xy-planc 
is calculated by (a) above and (16a) to be: 



dS= 11 yjy 2 + x 2 + 1 dx dy = f f y/r 2 + lrdrd6 = 2vr- : !-(r 2 +l) 3/2 

S J J R Jo Jo 
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